home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor2 / hamilton.src < prev    next >
Text File  |  1993-02-18  |  2KB  |  79 lines
  1. %%HP: T(3)A(R)F(.);
  2. DIR
  3.   SetUpHamiltonian
  4.    \<< 2 ^ '0.5*(px^2+py^2+pz^2)' + 'H' STO
  5.        H 'x' \.d COLCT NEG 'pxdot' STO
  6.        H 'y' \.d COLCT NEG 'pydot' STO
  7.        H 'z' \.d COLCT NEG 'pzdot' STO
  8.        H 'px' \.d COLCT 'xdot' STO
  9.        H 'py' \.d COLCT 'ydot' STO
  10.        H 'pz' \.d COLCT 'zdot' STO
  11.        'x' Lie4thOrder '\Gdx' STO
  12.        'y' Lie4thOrder '\Gdy' STO
  13.        'z' Lie4thOrder '\Gdz' STO
  14.        'px' Lie4thOrder '\Gdpx' STO
  15.        'py' Lie4thOrder '\Gdpy' STO
  16.        'pz' Lie4thOrder '\Gdpz' STO
  17.    \>>
  18.   RunHamiltonian
  19.     \<< 6 PICK 5 PICK R\->C
  20.       \-> U
  21.       \<< \Gd \Gd * 2 / DUP '\Gd2' STO
  22.           \Gd * 3 / DUP '\Gd3' STO
  23.           \Gd * 4 / '\Gd4' STO
  24.           { # 0d # 0d } PVIEW
  25.         DO
  26.           HamiltonStep
  27.           6 PICK 5 PICK R\->C
  28.           DUP U LINE 'U' STO
  29.         UNTIL KEY
  30.         END DROP {x y z px py pz} PURGE
  31.       \>>
  32.     \>>
  33.   HamiltonStep
  34.     \<< 'pz' STO 'py' STO 'px' STO
  35.         'z' STO 'y' STO 'x' STO
  36.         \Gdx \->NUM x +
  37.         \Gdy \->NUM y +
  38.         \Gdz \->NUM z +
  39.         \Gdpx \->NUM px +
  40.         \Gdpy \->NUM py +
  41.         \Gdpz \->NUM pz +
  42.     \>>
  43.   Lie4thOrder
  44.     \<< LieDer COLCT DUP '\Gd' *
  45.         SWAP LieDer COLCT DUP '\Gd2' *
  46.         ROT SWAP + SWAP
  47.         LieDer COLCT DUP '\Gd3' * ROT SWAP
  48.         + SWAP LieDer COLCT '\Gd4' * +
  49.     \>>
  50.   LieDer
  51.     \<< \-> N
  52.       \<< N 'x' \.d xdot *
  53.           N 'y' \.d ydot * +
  54.           N 'z' \.d zdot * +
  55.           N 'px' \.d pxdot * +
  56.           N 'py' \.d pydot * +
  57.           N 'pz' \.d pzdot * +
  58.       \>>
  59.     \>>
  60.   H 'px^2+py^2+pz^2+(x^2+y^2-z)^2'
  61.   PPAR {(-6.5,-3.1) (6.5,3.2) X 0 (0,0) FUNCTION Y }
  62.   xdot 0
  63.   ydot 0
  64.   zdot 0
  65.   pxdot 0
  66.   pydot 0
  67.   pzdot 0
  68.   \Gdx 0
  69.   \Gdy 0
  70.   \Gdz 0
  71.   \Gdpx 0
  72.   \Gdpy 0
  73.   \Gdpz 0
  74.   \Gd .1
  75.   \Gd2 .005
  76.   \Gd3 1.66666666667E-4
  77.   \Gd4 4.16666666668E-6
  78. END
  79.